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\ Abstract 



We study parametric integration of functions from the class 
C r ([0,l] dl+d2 ) to C([0, in the quantum model of computa- 
| tion. We analyze the convergence rate of parametric integration 

m 

O 



in this model and show that it is always faster than the opti- 
mal deterministic rate and in some cases faster than the rate of 



^ ■ optimal randomized classical algorithms. 

CO 

o 

1 Introduction 

Summation and integration are the most famous numerical problems that 
achieved a speedup in the quantum model of computation, compared to the 
optimal convergence rates of deterministic and randomized algorithms in the 
classical case. 

In this paper we study the problem of parametric integration, where the 
integral depends on a parameter. Therefore, the solution is now a function, 
so the problem carries features of both integration and approximation. 

We will consider the problem from the point of view of complexity theory 
and provide an analysis for the class of r-times continuously differentiable 
functions. For this class we determine the order of the minimal error (up to 
a logarithmic gap) by deriving matching upper and lower complexity bounds. 

In Section |21 we present the required notions from quantum information- 
based complexity theory, recall related previous results and formulate the 
main result. Section |3] is devoted to the proof of the upper bound. In 
Section |3] we prove the lower bound, and in the final Section El we give some 
comments on the results. 
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2 Preliminaries 



In this Section we formulate the problem which is investigated. Then we give 
the basic definitions of quantum information-based complexity theory, state 
some useful technical results and finally formulate the main result. 

2.1 Problem formulation 

Let D l = [0,l] dl and D 2 = [0, Vf 2 with d u d 2 > 0. On the domain D = 
D\ x D 2 we define the function class C r (D) for an integer r > 1 as the set 
of all functions f(s,t), for which all partial derivatives up to order r exist 
and are continuous. Let a be a multiindex, then the norm ||.|| r on C r (D) is 
defined by 

\\f\\ r :=m r \\f^\\ c{D) . 

\a\<r 

Let C(Di) be the space of continuous functions on D 1 with the supremum 
norm. We consider the solution operator 

S : (7(D) -> C(Di) 

f » (Sf)(s)= [ f(s,t)dt. (1) 

This means, we study parametric integration: Integrate the family of func- 
tions f(s,t) parametrized by s G D\ over t 6 D2. The limiting cases where 
either d± = (pure integration) or d 2 = (pure approximation) are formally 
included because they represent classical problems of numerical mathemat- 
ics. The aim of this paper is to study the intermediate cases where di 7^ 
and d 2 7^ 0. 

2.2 Quantum Setting 

We use the terminology developed by Heinrich in [3], which is a translation 
of information-based complexity (IBC) methods to the quantum model of 
computation. In order to be as selfcontained as possible, we summarize the 
quantum IBC notions needed in this paper. 

First, we briefly recall the standard notation of quantum computing. Let 
H\ be the two-dimensional complex Hilbert space C 2 and 

H m = H x ® • - • <g) H ! 
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be the Hilbertian tensor product of m copies of H 1 . We use the following 
notation, 

Z[0,N) := {0,...,N -1} 

for JV 6 N. Let C m = {\i) : i G Z[0,2 m )} be the set of unit basis vectors of 
H m , also called classical states or basis states, and let U(H m ) denote the set 
of unitary operators on H m . 

Let F{D, K) be the set of mappings / : D — > K. Now we introduce the 
notion of a quantum query. For F C T a quantum query is given by a tuple 

Q = (to, to', to", Z, r, (3), 

where to, to', to" G N, to' + to" < to, Z C Z[0, 2 m ') is a nonempty subset, and 



r : Z 
/3 : if 



Z[0,2 m ") 



are arbitrary mappings. The mapping r is the coding from basis states of 
one register of the quantum computer to the domain of / G F, whereas (3 
is the coding of the function values from the range of / to basis states of a 
second register of the quantum computer. Such a tuple Q defines a query 
mapping 



Q : F 
f 



U{H m ) 

Qf 



by 



Qf\m\y) 



\i)\x®p(f(r(i))))\y) if^z 



(2) 



\i)\x)\y) 



otherwise, 



where \i) G C m >, \x) G C m », \y) G C m -m'-m" (if m = to' + to", we drop the last 
component) and © means addition modulo the respective power of 2, here 
modulo 2 m . The total number of qubits needed for Q is m(Q) = to. 

Suppose we are given a mapping S : F — » G, where G is a normed space 
(in this context S is a general mapping). We want to approximate S(f) for 
/ G F with the help of a quantum computer. To do so, we formally define 
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the notion of a quantum algorithm. A quantum algorithm on F with no 
measurement is a tuple 

A=(Q,(U J )] =0 ), 

where Q is a quantum query on F, n G No and Uj G U(H m ) (j — 0, . . . , n), 
with m = m(Q). Given such an A and / G F we define Ay G U(H m ) by 

A f = U n QfU n - 1 ...U 1 QfU . 

By := n we denote the number of queries and by m(A) = m = m(Q) 

the number of qubits used by A. We also introduce the following notation. 
Let A f (x,y) for x,y G Z[0,2 m )be given by 

A f \y)= A ffay)\x)- 

xez[o,2 m ) 

Hence (Af(x,y)) XtV is the matrix of the transformation Af in the canonical 
basis C m . 

A quantum algorithm on F with output in G with fc measurements is a 
tuple 

where k G N, and A; (/ = 0, . . . , fc — 1) are quantum algorithms on F without 
measurement. We set mi = m(Ai). Then bo G Z[0, 2 m °) and for 1 < I < k — 1, 
bi is a function 

i-i 

h : JJZ[0,2 mi ) -> Z[0,2 mi ), 

i=0 

and (/? is a function with values in G, 

k-i 

up : JJZ[0,2 mi ) -> G. 

/=o 

The function </? combines the outputs of the algorithms to give a final 
result. The functions b\ determine the starting state of the next algorithm 
Ai depending on the results of the previous algorithms. 

We also say that A is a quantum algorithm with measurement (s), or just 
a quantum algorithm. 

Let Vq(G) be the set of all probability measures on G whose support is a 
finite set. The output of A on input / G F will be an element A(f) G Vq(G) 
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(we use the same symbol A for the mapping A : F — > P (G)). We define 
A(f) via a sequence of random variables (£;,/) zLo 1 (we assume that all random 
variables are defined over a fixed - suitably large - probability space (Q, E, P)). 
Let now / e F be fixed and let £jj be such that 

F(^ f = x) = \A OJ (x,bo)\ 2 

and, for 1 < I < k — 1, 

p (6,/ = z|fo,/ = x ,---,6-i,/ = ^-i) = \Aij(x,bi{x , . . . ,xi-x))\ 2 . 

This defines the distribution of (Cz,/)^ 1 uniquely Let us define for x e 
Z[0,2™°),...,x fe _ 1 eZ[0,2™*-i) 

PA,/(aJo, • • • ,Zk-i) = |^o,/(^o,&o)| 2 |^i,/(^i,&i(^o))| 2 --- 

• • • \A k -i,f(x k -i,bk-i(xo, . . . ,x k - 2 ))\ 2 - 

It follows that 

P(&,/ = X , • • • , ffc-i,/ = X fc _i) = PAj(x , • • • , X fc _i). 

Finally we define the output A on input / as 

A(/)=dist(y9(£ ,/,...,a-i,/)), 

the distribution of <p(£o,f, • • • , Ck-i,f)- 

The number n q (A) := n q (Ai) is called the number of queries used 

by A. This is the crucial quantity for our query complexity analysis. 

Now we define the error of a quantum algorithm A: Let < 9 < 1, and 
let C be a random variable with distribution A(f). Then the (probabilistic) 
quantum error of A for S on input / with failure parameter 9 is defined by 

e(S, A, f, 9) := inf {e | P(||5(/) - Clio > e) < 0}. 

Then we put 

e(S,A,F,9) :=supe(S,AJ,9) 
feF 

and 

el(S, F, 9) := inf{e(S, A, f, 9) \ n q {A) < n}. 
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We will consider these quantities for the fixed error probability 1/4 and set 

e(S, A, f) = e(S, A, f, 1/4), e(S, A, F) = e(S, A, F, 1/4), 

and we define the n-th minimal query error of the problem class F and the 
mapping S by 

el(S,F):=el(S,F,l/±). (3) 

This means that we will analyze the given cost. There is a close 

connection between e q n and the e-complexity of a problem, which is defined 
by 

comp^S 1 , F) := min{m | e q m (S,F) < e}. 

The two quantities satisfy the following relation: For all n G No, e > we 
have 

e q n (S, F)<e^ comp^S, F) := min{m | e q m (S, F) < £l } < n Vei > e. 
2.3 Tools from Quantum Complexity 

For our analysis of parametric integration in the quantum model we will need 
some statements from quantum IBC, which are now summarized: 

Lemma 2.1 Let F C T{D, R), I G N and let S k : F -> G (k = 0, . . . , I) 

be mappings, where G is a normed space. Define S : F — > G by S(f) = 
Yl l k=o S k(f)- Let$ o ,...,0 l > 0, n , . . . ,n z G N and put n = Ylk=o n k- Then 

i i 

el(S,F,J20k) <J2<( S k,F,9 k ). (4) 

fc=0 k=0 

This is a generalization of Lemma 2 from Heinrich [7] , which can be proved 
by the same technique by just replacing the absolute value with the norm on 
G. 

Lemma 2.2 Let f(D, K) and ^ F C T{D, K). Let V : F -> F 
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be of the following form: there exist k, m* G N and mappings 

Vj : D^D (j = 0,...,«-l) 
(3: K-^Z[0,2 m ") 

p: D x Z[0,2 m *) K -»• X, 

swc/i t/iat /or f £ F and s & D 

(T(f))(s) = p(s, /3 o / o Vo (s), . . . , (3 o f o Vk ^(s)). 

Given a quantum algorithm A from F to G, there is a quantum algorithm A 
from F to G with 

n q (A) = 2Kn q (A) 

and for all f G F 

A(f) = A(T(f)). 

Consequently, if S : F — > G is any mapping and S = S o T , then for each 
n G N 

el Kn (S,F)<eUS,F). (5) 

The proof of Lemma 12.21 can be found in Heinrich [7] . 

We finally state some calculation rules for the query error: 

Lemma 2.3 Let S,T : F — ► G be mappings, n G N and e^(S,F) be finite. 
Then it holds: 

(i) 

e n q (T,F) < e n q (S,F) +sup\\T(f) - S(f)\\. (6) 

feF 

(ii) For A G R it holds 

e n q (\S J F) = \\\e n q (S,F). (7) 

(Hi) If K = R and S is a linear operator from F(D,K) to G, then for all 
AeRwe have 

e n q (S,XF) = \X\e n q (S,F). (8) 



The proof of this Lemma can be found in Heinrich 0j. 

We now cite a method how to increase the success probability of a quan- 
tum algorithm. Let M G N and ipo '■ ^ M — » M be the median of M numbers. 
For a quantum algorithm A we define ipo(A M ) := ?po(A, . . . , A) to be the 
median of the results of M repetitions of A. 

Lemma 2.4 Let T : F —>■ R be a mapping and A a quantum algorithm. 
Then 

e (T, MA M ), f, e" A//s ) < e (t, A, f, ±) . (9) 
A proof of this Lemma can be found in Heinrich 0] . 

2.4 Main Results 

First we recall the known results for the special cases where either d\ = or 
d 2 = 0. 

We use the asymptotic notation a n =4 b n for sequences of nonnegative real 
numbers a n and b n , which means that there exist some constant c > and 
some n G N such that a n < cb n for all n > n . If a n =4 b n and b n =4 a n 
then we write a n x b n . We often use the same symbol c for possibly different 
constants. In particular, the needed constants may depend on d ly d 2 and r. 

Let 13(G) denote the unit ball of a normed space G, i.e. 

13(G) :={geG\ \\g\\ G <l}- 

Then we have the following two theorems that are important for our analysis. 
The first considers quantum integration and was proved by Novak |lUj . 

Theorem 2.5 Let S mt : C r (D2) — > K be the integration operator. Then 

e q n (S int , B (C r (D 2 ))) x rr r ^-\ (10) 

Let A int (-, n) be a sequence of quantum algorithms which is of optimal order, 
that is 

e(S in \A int (.,n),B(C r (D 2 ))) < cn r l d ^\ 
The second theorem is concerned with approximation. 
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Theorem 2.6 Let S appr denote the function approximation problem, that is 
the embedding operator from C r {D{) to C(Di). Then 



e7(^ appr ,i3(C r ( J D 1 )) x eZ(S appr ,B(C r (D 1 )) 

x e q n (S appr ,B(C r (D 1 ))^n- r/dl . (11) 

Here e det and e™ c are the minimal deterministic and Monte Carlo error. A 
detailed definition can be found in [S]. The proof of the rate for the quantum 
case is due to Heinrich 

We state another important result, which is in fact the key to the inte- 
gration result mentioned above. Let 

LZ ■= {</:{0,...,iV-l}-R}, 

with the norm ||<7||£jv = maxj |g(«)|. Then we get the optimal rate for quan- 
tum summation, where the upper bound is from Brassard, H0yer, Mosca and 
Tapp P and the lower bound from Nayak and Wu jH], with the extension 
that is used for our purpose coming from Heinrich jl]: 

Theorem 2.7 Let '■ L^ — > R be given by S^g = i J^iLo 1 9^)- Then for 
n < N there are constants C\, C2 not depending on N such that 

cm- 1 <el(S N ,B(L^))<c 2 n- 1 . (12) 

Note that the convergence rate does not depend on the number of summands, 
so we can choose iV to be large enough to satisfy any needed precision for 
the approximation of an integral by such a weighted sum. We denote by 
A sum {-,n, N) a sequence of quantum algorithms with this convergence rate, 
meaning that 

e(S N ,A sum (;n,N),B(LZ)) <c 2 n~ l . 

Finally we consider the rates of parametric integration. From now on 
let F denote the unit ball of C r (Di x _D 2 )- In order to enable comparison, 
we also recall the results in the deterministic and Monte Carlo setting. A 
detailed analysis for the Monte Carlo setting can be found in Heinrich and 
Sindambiwe jHj. 



9 



Theorem 2.8 The minimal errors of the operator S from (QJ) behave in the 
following way: 

-r/(d!+d 2 ) 



eT(S,F) 



\S,F) 



n 



n -(r+d 2 /2)/(di+d2)(l g n \l/2^ 

r / dl (logn) r / dl , 



and 



n 



-1/2, 



(log n) 1/2 e™(S,F) 4 n- l ''\\ogn 



1/2, 



\3/2 



ifr >d 1 /2 
ifr<d 1 /2 

ifr = d 1 /2. 



(13) 



(14) 



The following theorem settles the quantum query complexity of parametric 
integration and answers the question when and how much better quantum 
algorithms are (as compared to deterministic and Monte Carlo methods). 
The comparison is discussed in detail in Section El 

Theorem 2.9 For (ii ^ 0, ^ ^ and r > the operator S from (QJ) satisfies 



r + d 2 



r+d 2 



r+d 2 



n 



n 



^1+^ 4el(S,F)4 n d i+ d 2{\ogn) d ^- 



-r/di 



ifr>di (15) 
ifr<d l (16) 



and 



n 



-i 



4eUS,F)^n~ L (\ognY, if r = d x . 



(17) 



Note that in asymptotic statements we leave the logarithm unspecified, 
whereas in cases in which the basis is essential we write, e.g., log 2 n or Inn 
to indicate base 2 or the natural logarithm. The proof of the above theorem 
consists of two parts: First we prove an upper bound for the query error 
e^(S,F). Second, we prove a lower bound for e^(S,F) which has the same 
order as the upper bound, up to a logarithmic gap. 



3 Upper Bound 

Now we assume d\ ^ 0, d 2 ^ and r > 0. Let k G No be fixed, and let 
denote the partition of D\ into cubes of sidelength 2~ k with disjoint interior. 
Let 

A« := {r- 1 2- k {i 1 ,...,i dl ):0<i 1 ,...,i dl <r2 h } 
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be the equidistant mesh of sidelength r 1 2 fc on Dj. Define 

n 1 , k :=\A k 1) \ = {r2 k + l) d \ (18) 

Let 

P, (1) : C r (D 1 ) -> C(D0 

be the <ii-dimensional composite Lagrange interpolation of degree r on A^\ 
This means, on each cube Q G nj^ the function P^ f is the <ii-dimensional 
tensor product Lagrange interpolation over the nodes Q n A^. Note that 
the resulting function is an element of C(Di). 

For / fixed PjPf is uniquely defined by {/(s) : s G A^}. Therefore the 
operator P^ will also be interpreted as defined on L^A^), the space of 
real valued functions on A^, equipped with the maximum norm. 

Finally we also consider the operator P^ as acting in the space C(D 1 x 
D 2 ), meaning that we interpolate with respect to the first component only, 
leaving the other one fixed. In this case P^ is defined by (P^ f)(s,t) : = 

(n (1 V(-, *))(*)• 
Let 

that is .Di = Uj=o Qkj anc ^ the Qfc? are cu bes of sidelength 2 with disjoint 
interior. Let s k j be the point in Q k j with the smallest Euklidean norm. We 
define the restriction operator R kj : ^(D^M) — > ^"(Z?!, M.) by 



(R kj g)(s) = < 



5 (2 k (s - s fc j)) , if s e Qfej 

otherwise. 



(19) 



Let v = (r + l) dl and let 

i = h(r + if 1 ' 1 + i 2 {r + l) dl - 2 + ... + i dl _i(r + 1) + i dl 

for i = 0, . . . , v — 1 be the representation of i in base r + 1. Let now 0j (i = 
0, . . . , v — 1) be the tensor product Lagrange base polynomials of degree r on 
Di for the grid Aq 1 " 1 , meaning that <f>i(s) = 1 at the point s = r~ 1 (ii, . . . , i^) 
and (f>i(s) = for all other points in A^ . 
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Since we want to take advantage of the fast convergence of quantum 
summation, we have to find a function whose integral can be approximated 
by quantum summation. This function needs a small supremum norm to 
give quantum summation its full impact. For each level k > 1 and a fixed 
gridpoint s G \ Aj^ we define the detail function f kjS G C r (D 2 ) as the 

difference function between / and its approximation P^f, both functions 
considered for this fixed s. 

Let us consider the structure of the detail function in dependence of / and 
the tensor product Lagrange base polynomials. For fixed s G Aj^^A^ there 

is a cube Qk-i,j( s ) £ with s G Qk-i,j(s) (if there are several possibilities, 
choose the one with the smallest index j(s)). The detail function has the 
form 

v-1 

f k , s (t) = /(M) - (P&/XM) = f( s ,t) - Y^Rk-uwMWfat) , (20) 

i=0 

where 

Sl = s k - ltj{s) + ^-^(h, . . . ,i dl ). (21) 

The following Lemma shows that the detail function has a bounded ||.|| r - 
norm. 

Lemma 3.1 There is a constant c > such that for any function f G 
B(C r (D)), any integer k > 1 and any s G A^ \ Aj^ we have 

\\fk,s\\r<C (22) 

Proof: 

The functions Rkj<t>i (i — 0, . . . ,v — 1) are the tensor product Lagrange base 
polynomials on Q k j for the grid Qkj H and 

||0i||c(Di) = \\Rkj(t>i\\c{Q kJ ) (23) 
for k G N and % — 0, . . . , v — 1. Since sup^ ||0i||c(Di) < c' we get 

v-1 

\\fk,s\\r < \\f(s,t)\\ r + \(Rk-lj(8)<l>i)(s)\ \\f(Si,t)\\ r < 1 + Vc' . 

i=0 

We now choose c = 1 + v c' and the statement follows. □ 
We will need the following 
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Lemma 3.2 Let k > be an integer. Then the operator P fc , considered as 
acting from L^' k to C(Di) is bounded by a constant which does not depend 
on k. 

Proof: 

Let z G Lc£' k with \\z\\ T n i,k < 1. With the notation from above we infer that 

"J-'oo 

for j such that t G Qkj we have 

u-l 
i=0 

and therefore 

v-l 

||-Pfc (1) «l|t7(I3i) < Yl \ Z ^ \ ' H^i^llcCe^) <VC' < C, 
i=0 

where c is the constant from the proof of Lemma 13.11 which is independent 
of k. □ 



Let us now state the parameters that are needed for the proof of the 
upper bound. We use a multilevel approach developed by Heinrich (see 0) 
which was also used to obtain the optimal Monte Carlo rates for parametric 
integration in [S]. For x G R the notation \x\ means the smallest integer 
greater than or equal to, and [^J the greatest one smaller or equal to x. For 
neNwe set 

-(log 2 n) + l 



d 1 + d 2 

The starting level m is defined by 



(24) 



m 



and the final level I by 

r(l + d 2 /r)ml, 



rn, if r > di 
otherwise, 

if r > o?i 



(25) 



/ := < 



(26) 



[(l + ^/^i) 7 ^! — P otherwise, where p := [(log 2 m)/d\\ . 



13 



We use n it k points for the interpolation on level k, and we recall that n^ k = 
(r2 k + 1)*. Let 

M k := [8(A; + 3)ln2 + 81nn lifc l , (27) 
then we define the query number for quantum summation as 



n 2 ,k ■-- 



^2d 2 m-±(r+di)(k-m)^ if r > d\, 

[M^ 1 2 (dl+rf2)m ~ dl M( d i- r )( i - fe )] otherwise. 



(28) 



The number of summands for quantum summation in level k (k — m, . . . , I) 
is defined as 

N k := 2 rkd *n d 2 %. (29) 

Let us shortly describe the main idea of the proof: In the starting level we 
approximate those integrals directly, which correspond to parameters s on 
the roughest grid, with the finally needed accuracy. On the finer levels we 
do the same for the detail functions. Then we interpolate the computed 
approximations and add them up to get our approximation to the solution 
function. 

Now we prepare the discretization of the functions that will be used on 
the quantum computer. To do this, we need a mapping from our function 
class F to L^. For k > m we choose a number m* of qubits so large that 



and 

Then we define 
(3:R -> 



2 m*/2-l > 1 



2 -rn*/2 < 2- rk n~l. 



(30) 



{0,...,2 m *-l} 




if z < -2 m */2-i 



Z I— > < 



\2 m */ 2 ( z + 2 m */ 2 - 1 )J if - 2™*/ 2 - 1 < z < 2 m */ 2 - 1 



2 m - 1 



if z > 2 m */2-i. 
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Furthermore we define 

7 : {0,...,2 m * - 1} -> R 

y 2 m * /2 y -2 m * /2 -\ 

On the starting level m we only have to approximate the integral of / for 
fixed s G A~', so in this case we just discretize the function /(s, ■), which 
means that we have a function n : {0, . . . , Nfh — 1} — > x _D 2 which is 
defined by ^(j) = (s,tj), where the points i,- G -D2 are node points needed 
for quantum summation, they will be specified below. Thus for the starting 
level we get 

(W)(i)=7((/3o/o^0')))- (31) 

Let now k > in be fixed. As already indicated, we will approximate the 
integral of the detail function for fixed s G Aj^Ajj. 1 ^ by quantum summation. 
For a fixed summation number from (J29j) we define 

iji : {0, . . . , N k - 1} -+ Di x D 2 (i = 0, . . . , v) 

by 

Viti) ■= ( s i>*i) , 

where s v = s and the points S{ (i — 0, . . . , v — 1) are the points from (J5TJ). 
The points tj G -D2 are again node points needed for quantum summation. 
Finally we define 

p : {0,...2 m * - 1} U+1 -> R 

by 

u-l 

p(2/o, •..,!/«):= 7(?/«) - ^(-Rjfc-ij-W^OWT^i)- ( 32 ) 

i=0 

From these mappings we get the operator : F — > by 

r M / := p (09 o / o ^ =0 ) . 

This means that 

Now we are ready to compute the query error of S. By Lemma 12.31 we 
can decompose the query error into 

el(S, F) < sup \\Sf - P®Sf\\ + el(P^S, F). (33) 

feF 
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So the error splits into a deterministic and a quantum part. Classical polyno- 
mial approximation gives for g e C r (Di) and k G N , (see e.g. [2J, Chapter 
3.1) 

\\g-Pi 1) g\\c(D l) <c2- r \ (34) 
so for the deterministic part in ()33|) we get 

\\Sf-P l {1) Sf\\ c{Dl) <c2- rl . (35) 

Next we consider the quantum part of Let 

J = Jife*" 1 + J2b d2 ~ 2 + ...+ 3d2 -ib + j d2 

for j — 0, . . . , Nfe — 1, where b = 2 rk n 2t k- Let the node points for the quantum 
summation be defined as 



U 



3i 32 3d 



For k > in we define the operators Jk )S '■ F — > R by 

4, J := ^ ( 36 ) 
fe i=o 

which is the rectangle rule with points for ff. s . Next we define, also for 
k > m, operators Uk )S '■ F — ► R by 

U k J = [ h, s (t) dt. (37) 

Since the accuracy of the rectangle rule with N k points in dimension d 2 is of 
the order N^ 1 ^ 2 for functions with bounded first derivatives, we get by (}2"2"j) 
and pUj) 

\U k M)-JkM)\<c2- rk n~l. (38) 
By definition of the discretization operator )S we get for \z\ < 1 

7(^))<^<7(^)) + 2" m * /2 , 
and by (|3Uj) this implies that 

|r M (/)(j)-A.. (*i)l 

< -7W(s,^)))| + ^1^-1^(5)1 (/(^^-TW^.i;)))! 

i=0 

< c2-^n-i. (39) 
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From the discretization accuracy of )S and (|34j) we also infer that 

l|r*,.(/)|| L ^ < \\(f h ,.(tj))\\ L »k + \\(fk, s (tj)) - r M /ll^ 

< c2 ^( fe - 1 ) + cri^T Tk < Cl 2~ rk , (40) 

which implies 

r fe)S (F) C Cl 2-' k B (L^) . (41) 
From it also follows that 



\S Nk T k J ~ J k , s f\ < ^ E K r M/)(j) - < c2- rk n^ k . (42) 

k 3=0 

Now we calculate the error of the integration of the fk, s on \ A^. 
We get with Lemma ESI (021), Theorem E3 Lemma O and (jUJ) 

< SUp \U kta {f) - J*, a (/)| + 4(v+l)n 2 , k ( J k,s, F) 

< sup | £/*,,(/) - Jfc,»(/)| + sup |SW fc rfc,./ - Jfc jS /| + e^ w+1 , (^7v fe r fcjS , F) 

< c2- r %' + el 2 k (S Nk , Cl 2- rk B(L% ) ) 

< c2^ k n-l + Cl 2- rk el 2k {S Nk MLl k )) 

< c2~ rk n^ k + c2- rk n^ k 

< c2- rk n~l. (43) 

With the help of this result we can now investigate the error of the operator 
P^S. Since 

P?>S = PUS + £ {PP-P&)S, (44) 

k=m+l 

we investigate the error of the operator 

(P^-P^S-.F^Cm. 

We define 

e k :=2- {k+3 \k = m,...J. (45) 
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Then we set 

h m := Mf h nx^n 2 ,ffi, (46) 
and for k = rh + 1, . . . ,1 we set 

h k := M k (n 1>k - ni tk -i)2(v + l)n 2}k . (47) 

Let A ks be a quantum algorithm that computes an approximation to 
U ktS on F with the rate from (J43j) and let ( ktS be a random variable with 
distribution A kjS . We define a random variable £ k with values in L^' fc as 
follows: For s G \ A.^_i we ^ 6c ( s ) be the median of M k independent 
copies of (k,s, that is, we repeat A ktS M k times. For s G k k _ x we set £k(s) := 
0. Since (U k>s f)(t) = for s G A^, by this choice we establish an error of 
zero in these points. We have 

" P&){Sf) = P k W ((/ - Pffx)^/)) = PPiUH,)^, (48) 

where on the right hand side pj^ is considered as acting on L^' fe . This 
means that because we can interchange interpolation with respect to the 
first component and integration with respect to the second component, we 
indeed compute an approximation to (P^ — P k ^i)(Sf) by (i^ — P k ^i)^ k - 
By Lemma and (tOl . 

P(|t/ M (/) " M\ > c2' rk n^) < e~ M *l\ (49) 

Consequently, 

n\U k , s (f) - < c2~ rk n 2 l V a G A«) 

> 1 - n ljfc e- Mfc/8 > 1 - 2~ (fe+3) = 1 - 9 k (50) 

by (J23)- From Lemma E!2J (|4T|) . (}4*%j) and (j^U)) we obtain for the query error 
of the operator (P k ^ — Pf^AS, 

4 (Vi 1} - f £i)< s ', F ' **) ^ c2 ~ rfc ^- (si) 

We use Lemma ETTl to calculate the error of P^S. From © we get Y! k=m &k < 
1/4, hence with Lemma 12. 11 Lemma [3. 21 (J44)) and with 

l 

n := n fc , (52) 

fc=rrt 
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we get 

el(P^S,F,l/4) < e^P^F,^) 

+ E 4M 1] - p k-^ F ^^ ( 53 ) 



k=m+l 

Now we consider the different cases: 



1. r < d\. For the error on the starting level m = we can make direct use 
of Theorem 12.71 and, in this case using the operator from (|31|). we get 

by Theorem 12.71 Lemma 13.21 and a similar probability argument as above. 
Now with (fo^j) and (foTj) we get 

4 (p/%,F,1/4) 

< c(n 2 , )- 1 + ^c2-^n 2 J 

k=l 

I 

< c(log 2 n)^2-( dl 

fc=0 

i 

< c(\og 2 n)2-^ +d2)m J2 2 ^ dl 



i 

l,k 
k=l 

I 

-(di+d 2 )m+(iifc+5(cii-r)(J-fc)-rfc 



k=0 

I 

_ -r)(l+k) 
k=0 

From this we get with the help of the geometric sum formula 
eg (p/ 1} 5,F,1/4) < c(log 2 n)2- (dl+d2)m+(dl - r)z 

and with 2 dl '- rf2m x 2 dim - lo ^ m we arrive at 

e? (p/^F, 1/4) < c2~ H < cn- r/dl (\og 2 n) r/dl . 

For the deterministic part of the error we get by the choice of I 

\\Sf - P^SfW < c2- rl < cn- r l dl i\og 2 n) r / d \ 
which by (I33j) gives the desired rate for r < d\. 
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2. r > d\. To calculate the error on the starting level m = m we use Theorem 
12.51 and (J28j) and again with the probability argument from above we get 

eljpWS, F, 6 m ) < c (n^n'*- 1 < ^-C^)/^). (54) 

With flE3l and JSU we get 

< m -('"+ rf 2)/( rf l+ rf 2) + ^ C2 



I 

-rk -1 



fe=m+l 

(r+c(2)/(rfi+rf2) I „ o-rfc9-d2m+i(r-+di)(fe-m) 



< cn -( r + rf 2)/( rf i+ c( 2) + c ^ 2~ rfc 2 

fc=m+l 

Z 

< C ^~( r + c(2 )/( rfl + c( 2) _|_ c 2-( r + rf 2)"l ^ 2~5( r_d l)( fc_m ) 

fc=m+l 

I 



fc=m+l 

For r > di the sum is bounded by a constant, and for r = d\ the sum gives 
an additional factor of logn. By the choice of I we get 

- P^SfW < c2- rl < cn r { ' r+da)/( - dl+da \ 

and with (J33)) we arrive at 

e q h (S,F) < cn -(-+ d 2)/(rfi+^) (55) 

for r > di and 

el(S,F) < en' 1 logn (56) 

for r = d\. 

Finally we estimate the number of queries n that are needed to obtain 
the desired precision. Since the total number of queries is 

i 

h = M fh n lt f h n 2: fh + 2(v + 1) y~] M k (n hk - rii )fc _i)n 2 , fc , 

fc=m+l 



20 



we get for r < d\ 

i 

h = M n lfi n 2fi + 2(v + l)^2M k (n ltk -n lik „ 1 )n 2! k 

fc=l 

i 

< M k n ltk n 2 ,k 

fc=0 

i 

< c^2 dlfc 2 (dl " 



fc=0 

i 

+d 2 )m-dik-±(d 1 -r)(l-k) 



_ 2 

k=0 



fc=0 



■r)(i-fc) 



and for r > d\ we get 



n = M m n 1<m n 2<m + 2(v + 1) ^ M k (n 1>k - n 1)k -i)n 2 ,k 

k=m+l 

I 

< c ^2 M kn 1<k n 2jk 

k=m 
I 

< c^2(k + lnn hk ){r2 k + l) dl 2 d2m -^ r+dl ^ k -' m) 

k=m 

I 

-)d2rn—^(r+d\)(k—m) 



< clogn^2 dlfc 2 c 

k=m 

I 

< clogn 2 (dl+d2)m 2-^ {r ~ d ^ k - m) x n log ri, 



k=m 



since g^A; + d 2 m — i(r + di){k — m) = (d\ + d 2 )m — |(r — <ii)(fc — m). This 
means for r < d\ the cost is of order n and for r > d\ it is of order nlogn, 
such that a rescaling of n leads to the proposed rate. In the case r = d\ the 
cost is C(n(logn) 2 ), so together with we get the additional log-factor in 
the convergence rate. Now the upper bound of Theorem 12.91 is proved. 



Note that the proof of the upper bound was carried out in terms of query 
errors and can easily be translated into an explicit quantum algorithm for 
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parametric integration. The algorithm has the following form and uses the 
sequences of optimal algorithms A sum (-,n, N) and A int (-,n) for quantum 
summation and quantum integration. For a given n we recall the needed 

-(log 2 n) + 1 , starting level fh :— m if 



parameters, which are m = 

r > di and zero otherwise, final level / := [(1 + d 2 /r)m] if r > d\ and 
(1 + d 2 /di)m~\ — p otherwise, where p := [(\og 2 m)/di\. We have n ljk = 
(r2 k + l) d \ M k := \S(k + 3) In 2 + 8 lnn ljfc ] , and n 2 , k ■= \2^ m ~^ r+d ^ k - m '>], 
if r > d x and ^ M k ~ 1 2 ( - dl+d ^ m ~ dlk ~^ dl ~ r ^ l ~ k ^'\ otherwise. Finally we have 
N k := 2 rkd2 n d2 k . Now the algorithm A parint (f,n) is the following: 

1. Starting level fh: For all s G do: 

(a) If r > di, compute times A^ iS := A int (/(s, •), n 2 ,fh) and let 
Cm(s) De the median of these results 

(b) If r < d ± , compute times A^ s := A sum (T '^J ',n 2 ,m, Nf- n ) and 
let £m(s) be the median of these Mf^ results 

2. Finer levels: For k — fh + 1, . . . , I do: 

• For all s G \ A^i do: 

M fc times compute A Ks := c x 2~ rk A sum (c^ 1 2 rk T ktS f,n 2jk ,N k ) and 
let ^fc(s) be the median of these M fc results 

• For s G A^ do: &(s) := 

3. Final approximation: A parint (/, n) := ^~ 1) Cm+ELm+i( P fc 1)- - P fc-i)0 

In step 2 the function T k;S f G L^ fc is scaled by c[ 1 2 rk and the result then 
rescaled to make sure that the algorithm A sum is applied to a function with 
L^-norm smaller or equal to one. 



4 Lower Bound 

In this Section we first cite a general result for lower bounds on the quantity 
e%(S, F) and then we apply this result to the case of parametric integration. 
Let D and K be nonempty sets, let L G N and let to each 

u = (iio,...,UL-i) e {0,1} L 
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an f u G J-(D, K) be assigned such that the following is satisfied: 
Condition (I): There are functions go, gi G J-(D,K) and a decomposition 



D = \Jt=o A with A n D v = (/ ^ /') such that for t G D 



/«(*) 



5-0 (t) if u t = 
gri(t) if it; = 1. 

Next we define the function p(L, I, I 1 ) for L G N, < / ^ V < L by 



Note that j(L—j) is minimized iff \L/2—j\ is maximized. For u G {0, 1} L 
we set | it | = Yla=Q u i- Then we have the following 

Lemma 4.1 There is a constant cq > suc/i t/iat t/je following holds: Let 
D, K be nonempty sets, let F C J~{D, K) be a set of functions, G a normed 
space, S : F — > G a function, and L G N. Suppose (f u ) u e{o,i} L ^ ^ A i^) 
a system of functions satisfying condition (I). Let finally < / ^ I' < L and 
assume that 

f u G F whenever \u\ G {/,/'}• 

Then 

el(S,F) > ~rmn{\\S(f u )-S(f u ,)\\ : | u | = I, \u'\ = I'} 
for all n with 

n < c p(L, I, I'). 

A proof of Lemma 14.11 can be found in Heinrich ||. With the help of this 
Lemma we can now prove the lower bound for parametric integration: 

Let Co be the constant from Lemma [4.1 [ and let d = d± + d 2 . For n G N 
we choose an even number m G N such that 



L = m d > \J^ 2 n 2 + 4. (57) 

Let 

i = j x m d ~ l + j 2 m d ~ 2 + . . . + j d -im + j d 
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for i — 0, . . . , L — 1. Let 



ji + 1 



m m 



L-l 



, D = [J A 



j=0 



Define the functions ^ e C°°(.D) by 



. . . , S dl ,t dl+ i, . . . ,t dl+ d 2 ) — i)f\si, ■ ■ ■ , SdJ • V^^di+l, • • • , idi+d 2 ) 

dl d\+d,2 

= \\ri(msi-ji)- \\ rj(mti-ji) 

1=1 l=di+l 

with 



',(2)/ 



7}{x) = < 



e ^i^) ,0 < x < 1 
, otherwise. 



For k = Y?i=i d2 h ^ holds 

ftc kl ftc k2 f)f kd l+ d 2 

us 1 us 2 . . . oi di+d2 



1=1 l l=d\+l I 

di di+d2 

Y[v (kl) (msi- ji)m kl ■ |[ rf kl \mti- ji)m kl 

1=1 l=di+l 
di di+d 2 

m k Y[rj {kl \msi -ji) ■ Y[ r] {kl \mti- 

1=1 l=di+l 



From this we get 



C r (D) — max SU P 
Ztll d ' 2 h<r (s,t)eD 



d k iJi(s,t) 



ds kl ds k2 f)f kd i +d z 

Ub l Uh 2 ... OL dl+d2 
d\ di+d-2 

= max m k \\ sup |77 (fc ,)(s,)| ]) sup \rj (kl) (ti)\ 
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where 



di 



7 



d-i+d.2 

max If sup \v(ki)(si)\- ] J sup 



Therefore, setting ^ := ^Vi> we nave /« := Hi=o u i^i e F for a11 
Wj G {0,1}, i = 0, ... , L — 1. Since the ipi have disjoint support, f u satis- 
fies condition (I). Let cq = f Q r](x)dx > 0. It is easy to show that 

f a i 
\ ipi(si, • • • , td 1 +d 2 ) dsi . . . dtd 1 +d 2 = — -. 

J D m 

Let now I := L/2 - 1 and /':=/ + 1 = L/2. Then with (jSTj) it follows that 



p{L,l,l') > min y/j(L-j) 
j={i,i'} 



L 



- 1 



L 2 
T 



1 > c n. 



Since 



mm 



1 a, 



<h 



c(Di) r )m r m d2 



i>o (t dl+ x, . . . ,t dl +d 2 ) 



C(D!) 



d 1 +d 2 



= C n \ d i+ d 2^ ^+d 2 j ^ 

from Lemma l4~T1 we get a general lower bound for parametric integration and 
the lower bound of Theorem 12.91 for r > d\. 

We get a lower bound of n~ r / dl for approximation by simply choosing 
g?2 = 0, and since parametric integration can never have a better convergence 
rate than approximation, this gives us the lower bound of Theorem 12.91 for 
the case r < d%. 



5 Comments 

We have solved the problem of the quantum complexity of parametric inte- 
gration for the class C r (D) by providing upper and lower bounds, where the 
rates match up to a logarithmic factor. Now we compare our results to the 



25 



known results for deterministic and Monte Carlo methods stated in Section 
12.41 We again assume r > 1. Then the optimal quantum rates are always 
better than the optimal deterministic rate, except for the case c?2 = 0, where 
we have equal rates, and the following table provides a comparison of the 
deterministic, Monte Carlo and quantum case (without log- factors). 





e%\S,F) 


e™(S, F) 


e q n(S,F) 


r > d\ 


n ~r/(di+d 2 ) 


n -(r+d 2 /2)/(d 1 +d 2 ) 


n -(r+d 2 )/(di+d 2 ) 


di/2 < r < d x 


n ~r/(d 1 +d 2 ) 


n -{r+d 2 /2)/(d 1 +d 2 ) 


n -r/di 


r < di/2 


„— r / \d\+d 2 ) 


„—r/di 


n ~r/di 



We have to distinguish three different situations, depending on the relation 
of the problem parameters r and d\ . For r < d\ /2 the quantum rate provides 
an improvement over the deterministic rate, but it is as fast as the Monte 
Carlo rate. For this parameter constellation both the Monte Carlo and the 
quantum algorithm achieve the optimal rate of approximation which is the 
same in both cases. 

When we have di/2 < r < di, then the quantum rate is still the optimal 
rate of approximation, but the Monte Carlo rate is slower, which leads to the 
superiority of the quantum rate for this situation. 

For the case r > d\ we still have a better performance of the quantum 
algorithm as compared to Monte Carlo. 

Summarizing the discussion we can say that the quantum rate is always 
better than the deterministic rate, it is always at least as good as the Monte 
Carlo rate, and for r > d\/2 it is better than the optimal rate of Monte Carlo 
algorithms. 

Let us finally consider our problem in the bit model. Referring to jl] and 
[7j, we find that for our algorithm the number of qubits needed is (9(logn), 
the number of quantum gates is 0(n\ogn) and the number of measurements 
is C((logn) 3 log log n). 
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